Noise-induced dynamical phase transitions in long-range systems 
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In the thermodynamic limit, the time evolution of isolated long-range interacting systems is 
properly described by the Vlasov equation. This equation admits non-equilibrium dynamically 
stable stationary solutions characterized by a zero order parameter. We show that the presence of 
external noise sources, like for instance a heat bath, can induce at a specific time a dynamical phase 
transition marked by a non-zero order parameter. This transition corresponds to a restoring of the 
full ergodic properties of the system and may be used as a distinctive experimental signature of 
the existence of non-equilibrium Vlasov-stable states. In particular, we evidence for the first time 
a regime characterized by an order parameter pulse. Our analytical results are corroborated by 
numerical simulations of a paradigmatic long-range model. 



The study of the nonequilibrium behavior of long-range 
interacting systems recently triggered an intense research 
activity [l|. This is partly a consequence of the vari- 
ety of physical situations characterized by long-range in- 
teractions (including gravitational, plasma and nuclear 
physics, hydrodynamics, wave-matter interactions,...) 
[l|, which makes the topic amenable to interdisciplinary 
approaches and cross-fields fertilization. However, in 
most of the literature in this area the influence of an ex- 
ternal environment that may act as a thermal bath has 
not been taken into account. Especially for laboratory- 
scale systems, such an influence is expected to be rele- 
vant. Here we show that the existence of external noise 
sources can induce, on long-range systems characterized 
by different equilibrium phases, a dynamical phase tran- 
sition from a homogeneous Vlasov-stable phase to an 
inhomogeneous one. This dynamical transition corre- 
sponds to a restoring of crgodicity realized by the action 
of the noise. By extending the linear stability analysis 
of the Vlasov equation in order to take into account the 
influence of the noise on the system's dynamics, we iden- 
tify the time t c at which the dynamical phase transition 
occurs. At t c , in the thermodynamic limit, phase func- 
tions like the order parameter are not analytic. We thus 
provide a "dynamical phase diagram" in which time re- 
places the role of a thermodynamic parameter in ordinary 
equilibrium phase diagrams. For specific ranges of initial 
energies, the system undergoes a dynamical phase tran- 
sition even if the final equilibrium phase is homogeneous. 
This behavior corresponds to an order-parameter pulse 
occurring at t = t c . Our analytical results are corrobo- 
rated by numerical simulations of a paradigmatic model, 
which fully agree with the analytical predictions. 

A system is considered "long-range interacting" if the 
interparticle potential V decays at large distances r 
slower than l/r d , where d is the system's spatial di- 
mension. Under these circumstances, phase transitions 
have to be addressed within a proper definition of the 
thermodynamic limit, which is attained by rescaling 



the interaction strength with a function of the system's 
size in order to make the energy formally extensive [l[ 
(Kac's prescription) and by taking the number of par- 
ticles N —> oo. In this limit, the interparticle correla- 
tions (due to fmite-Af effects) are completely negligible, 
and the system's time-evolution is correctly described 
by the Vlasov equation [l|. We have recently demon- 
strated that a Hamiltonian thermal bath Q and a 
stochastic Langevin bath [|| produce the same equilib- 
rium and nonequilibrium effects when in contact with a 
long-range system, a conclusion not obvious a priori due 
to the character of the system's interactions. Assum- 
ing d = 1 for simplicity, the time evolution of the one- 
particle distribution function f(x,v,t) [normalized such 

that f_™ dx dv f(x, v, t) = N, Vt <G R] is thus con- 
veniently described by adding to the Vlasov equation a 
diffusion and a damping term. This gives the so-called 
mean-field Kramers' equation Q: 
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where ®{x,t) = /+~ dx' V {x - x 1 ) /+~ dv f(x',v,t), 
and D and 7 are respectively the diffusion and damp- 
ing coefficients. The former measures the strength of 
the noise perturbing the Vlasov dynamics, whereas the 
latter quantifies the dissipation to the external environ- 
ment. Below, we will assume that these two effects have 
the same physical origin, like when the system interacts 
with an heat bath. In particular, we suppose that the 
temperature T of the bath satisfies the Einstein rela- 
tion D — 7T (we work in natural adimcnsional units 
in which the Boltzmann constant and the mass of the 
particles are equal to one). One can however conceive 
more general situations in which the noise strength is 
not directly linked to the dissipation rate. An interesting 
recent example @ is the consideration of a noise term 
representing energy-preserving inter-particles fast colli- 
sion processes that are superposed to the long-range in- 
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teractions. We are interested in situations in which the 
potential V is such that there exist a critical temper- 
ature T c that separates a high-temperature equilibrium 
homogeneous phase from a low-temperature inhomoge- 
neous one. These phases can be described in terms of an 
order parameter, ip, which is zero at and above T c , and 
non-zero below T c . Furthermore, we assume that the 
interparticlc potential is even in the spatial coordinate: 
V(-x) = V(x). 

Let us consider a spatially homogeneous initial con- 
dition f(x, v, 0) = fo(v) that is Vlasov stable. For in- 
stance, this could be the result of a violent relaxation, as 
in Lynden-Bell's theory @. Hence, d^/dx = (<p = 0) 
and the system evolves under the sole effect of the inter- 
action with the thermal bath according to 
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Notice that for 7 = D = (microcanonical setup) , we re- 
cover the Vlasov stationarity condition df /dt = 0, which 
would last for infinite time in this thermodynamic limit. 
In contrast, for 7 ^= 0, the distribution will change with 
time due to the interaction with the bath. Equation ([2]) 
remains valid only as long as f(v, t) is dynamically stable. 
If, at time t = t c , it becomes Vlasov unstable, a dynami- 
cal phase transition occurs from the homogeneous ip = 
phase to an inhomogeneous ip 7^ phase [7( ■ The time t c 
can be identified by performing a linear (Landau) stabil- 
ity analysis [3, Q of the Vlasov equation for the f(v,t) 
which is solution of Eq. ^ at time t (in this analysis, 
f(v,t) is regarded as a function of v, the time t being 
fixed adiabatically). While a general f(v,t) can be stud- 
ied within the Landau formalism, in the following we will 
restrict ourselves, for the sake of simplicity, to distribu- 
tions with a single maximum at v = 0. In this case, one 
then shows 0, Q that t c is the minimum time t for which 
the equation 
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admits a solution for at least a k G M, k ^ 0. This ex- 
presses marginal stability. Here, V(k) £ R is the Fourier 
transform of the potential: V(k) = f^^dx e %kx V(x). 
In order to find t c , we can analytically solve Eq. ([2} for 
a given uniform initial condition f(x,v,0) = fo(v), plug 
the result in Eq. ([3]) , and invert the solution correspond- 
ing to the minimum time. 

The Green function of Eq. © is 

W(v, t\v , t = 0) = , (4) 

where a(t) = e"T*, and cr 2 (i) = T (l - e^'). Formally, 
the expression for f(v,t) is thus 
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dv o W(y,t\vo,t = 0)f o (v o ). (5) 



Since the time evolution can be expressed in adimensional 
units through the mapping t 1— > jt, below we implicitly 
adopt this convention. Our problem can be solved in 
Fourier space. First we take the Fourier transform in the 
velocities, [~](u>) = j+™ dv e lwv [•](«), of Eq.©, to get 



f(w,t) = f (a(t) w) e - 2 (*)- 2 /2. 



(0) 



Recalling the expressions for a(t) and <J 2 (t), it is clear 
that the short-time behavior of f(w,t) is dominated by 
the initial condition /o (w) (e.g., by the initial system's 
energy), while its long-time behavior is determined by the 
heat bath temperature. Inserting the Fourier representa- 
tion of f(v, t) in Eq. ([3]), we get that the transition time 
is equivalcntly identified by the zeroes of the function 



I(t,k) = 1 



V(k) 
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for some k ^ 0. 

The analysis outlined so far is valid under the follow- 
ing general assumptions: (i) One-dimensional system; (ii) 
Kac's thermodynamic limit; (iii) Symmetric potential; 
(iv) Spatially homogeneous distribution function f(v,t) 
with a single maximum at v = 0. We now specify the 
above general formulas to the Hamiltonian Mean Field 
(HMF) model, a long-range interacting system that has 
been widely used in the recent literature as a paradig- 
matic case of study [l|. The HMF model, which can 
be thought of as a set of N globally coupled rotators or 
XV-spins, is described by the Hamiltonian 
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where 9i € [0, 2tt[ is an angular coordinate, V(6) = 
(1 — cos9)/N is the (cosine) potential, and the order 
parameter p amounts to the system's magnetization: 

ip = ^2^ =1 {cosxi, sinxi) /N. Notice that, following the 
Kac's prescription, the interaction strength has been 
rescaled by a factor 1/N. In this way the total energy 
of the system, E = v i/^ + N(l - (p 2 )f2, becomes 

extensive, and we can consistently introduce the specific 
energy e = E/N. Once coupled with a Langevin ther- 
mal bath as in Eq. (TTJ), the system has also been named 
Brownian Mean Field (BMF) model 0]. 

To begin with, we consider uniform initial conditions: 

fo(v) = -^[&(v + v)-d(v-v)} (9) 

47T1> 

($ being the Heaviside step function), where v = 
V6(e - 1/2) and e is the initial specific energy of 
the system. In Fourier space we thus have fo(w) = 
Nsm{vw)/2TTvw, V(k) = tt (2S k ,o - 6 k ,i ~ /N, 
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FIG. 1: Stability analysis in the presence of external noise. 
Depending on the initial energy and on the temperature, an 
initially stable uniform state may stay uniform forever (a); 
becomes non- uniform at t — t c and remains non-uniform (b); 
becomes non-uniform at t = t c and then becomes uniform 
again (c). 



and we are led to examine the zeroes of [1 
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FIG. 2: (a) Magnetization pulse revealing a dynamical phase 
transition. Different curves are obtained for T = 0.6 > T c , 
N = 10 6 , averaged over 20 runs with initial conditions as 
in Eq. @. From top to bottom e = 0.60, 0.63, 0.65 < 
e c (0.6) = 0.662 . . .. (b) Finite-size analysis for T = 0.40 < T c 
and eo = 1.0. From top to bottom the lines refer to N = 
10 7 , 10 6 , 10 5 , 10 4 . 



I(t,±l) = l 
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If eo < e* = 7/12, 1(0, ±1) > 0, which means that at 
t = fo(v) is dynamically unstable [l], Since our 
analysis applies to states which are initially stable, in the 
following, we will restrict ourselves to cases with eo > e* . 
According to the values eo and T, there are three possible 
qualitative behaviors of I(t, ±1), as depicted in Fig. [T] 
In Fig. []Ji, I(t, ±1) < Vt > 0, and the system magne- 
tization ip is always zero (no dynamical phase transition 
occurs). Correspondingly, tc = +oo. This happens if T 
is above the critical value T c = 1/2 [J(+co, ±1) < 0] 
and if the initial energy eo is sufficiently high. An al- 
ternative is that I(t, ±1) crosses zero for a single value 
t = t c (Fig. [TJd). For t < t c <p is zero, whereas for t > t c 
the system becomes and remains magnetized. This oc- 
curs for T < T c . There is however a third possibility: 
I(t, ±1) may cross zero twice (Fig. [TJ;). This occurs for 
T > T c and if eo is below a critical value e c (T), which 
in general also depends on the shape of fo(v). The first 
zero corresponds to a dynamical phase transition occur- 
ring at t = t c , which drives the system to if > 0. Since 
T > T c , the system must eventually come back to zero 
magnetization in order to reach equilibrium. Under these 
conditions, one observes a "magnetization pulse" in the 
behavior of ip(t) (Fig [2k,). Notice that we cannot inter- 
pret the second zero of I(t, ±1) as the second transition 
time, since our stability analysis is only valid under the 
assumption of spatially homogeneous conditions, which 
breaks at t = t c - However, the system eventually goes 
back to ip — to equilibrate with the thermal bath. No- 
tice also that in this thermodynamic limit description, 
whenever the dynamic phase transition occurs, ip(t) is 



non- analytic at t — t c . This can be checked numerically 
with the BMF model, recalling that <p in a homogeneous 
state scales with the system size N as <p(N) ~ 1/yN 
(in the inhomogeneous phase, ip ~ 1). Rescaling by ^/N 
the curves obtained at different iV's, we indeed find that 
all the curves collapse for t < t c and develop the same 
non- regular behavior at t = t c (Fig. [5Jd). 

All the previous findings can be summarized in a dy- 
namical phase transition diagram representing t c vs eo, 
for different values of T (solid lines in Fig. [3j) . If 
T < T c , a dynamical phase transition always occurs 
and, correspondingly, < t c < +oo Veo > e*. For 
T > T c , the transition line ends at a critical value e c (T): 
if eo > e c (T), t c becomes infinite and the dynamical 
phase transition does not occur anymore. In the range 
e* < eo < e c (T), the dynamical evolution of tp corre- 
sponds in fact to a "pulse" starting at t = t c (see Fig. 
Ek)- Fig. [3] also reports the results of numerical simula- 
tions of the BMF model [H, which fully agree with the 
analytical curves. It is also possible to obtain two asymp- 
totic expressions for t c in a closed form [ll|: for T < T c 
and eo 3> 1, t c grows logarithmically as 
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For T>T C and eo — > e* + , t c goes to zero linearly as 
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We have also checked the validity of our conclusions for 
a class of (spatially uniform) initial distributions fo(v) 
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FIG. 3: Phase diagram of the dynamical phase transitions. 
At T < T c = 1/2, t c is always finite (plots refer to T = T c = 
0.5, T = 0.45, T = 0.4). Inset: At T > T c the curves end 
at a critical value e c (T) [plot refers to T = 0.55, e c (0.55) = 
0.584 . . .]. Above this energy, t c becomes infinite. Points and 
error bars refer to simulations with N = 10 7 , averaged over 
20 runs. 



wider than the water-bag, namely: 




x[&(v + v)-<&(v-v)]. (13) 

Here, 0<^<l/2wisa form factor with dimension of 
inverse velocity, and v depends now on both eo and k. 
For k = l/2v, Eq. (fl~3| reduces to Eq. ([9]), whereas 
for k = fo(v) becomes a bump- like distribution that 
touches the u-axis. Eq. (|13p can be understood as a way 
of modeling the action of a velocity-selector with milder 
cutoff effects than the water-bag, which are parametrized 
by the form- factor k. With respect to this class of ini- 
tial conditions, we have found that t c (n) is in general 
a monotonically decreasing function, meaning that for a 
given eo the water-bag initial conditions correspond to 
the most stable distribution of the class. However, all 



the qualitative features of Fig. [3] arc confirmed for any 
value of allowed k. 

The dynamical evolution of long-range interacting sys- 
tems is dominated by a global interaction which erases 
the interparticles correlations when N — > +oo. As 
a consequence, the system gets easily trapped in non- 
equilibrium states which correspond to the occupation 
of a small fraction of the available system's phase space 
0- We have shown that the action of a reservoir or of 
an external environment, renders these non-equilibrium 
states dynamically unstable, restoring the full (ergodic) 
occupation of the system's phase space. Correspondingly, 
the system manifests a dynamical phase transition, which 
can be detected by monitoring the time evolution of the 
order parameter. The transition time t c can be exactly 
identified by extending the Landau stability analysis to 
take into account the influence of the heat bath. The 
result is a dynamical phase diagram, in which time re- 
places the role ordinarily played by a thermodynamic 
parameter. The dynamical transition may occur even 
if the external temperature is above T c , and under these 
circumstances it is revealed by a peak in the time evo- 
lution of the order parameter. We have studied situa- 
tions in which the external noise satisfies a fluctuation- 
dissipation relation which enables the definition of the 
heat-bath temperature. However, it is also possible to ex- 
tend the Landau stability analysis to more general cases. 
This could be important to reproduce experimental con- 
ditions in which many sources of noise that do not nec- 
essarily satisfy a fluctuation-dissipation relation act on 
the system. The generality of our approach and the ro- 
bustness of our results under different initial conditions 
suggest that these dynamical phase transitions should be 
common and marked enough, to overcome experimental 
test. In particular, a pulse in the time evolution of the or- 
der parameter may be used as a signature of the existence 
of Vlasov stable out-of-equilibrium states in experiments 
with long-range interacting systems. 
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